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We study transport properties of an arbitrarily shaped ultraclean graphene sheet, adiabatically 
connected to leads, composed by the same material. If the localized interactions do not destroy 
chiral symmetry, we show that the conductance is quantized, since it is dominated by the quasi one- 
dimensional leads. As an example, we show that smooth structural deformations of the graphene 
plane do not modify the conductance quantization. 
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I. INTRODUCTION 

In recent years an explosive amount of attention has 
been paid to the novel material graphene^, a 2d carbon 
sheet. One of the important reasons for this great inter- 
est is the peculiar behavior of the associated transport 
properties. 

Essentially, graphene is a semi-metal with a gap over 
almost all the 2d Brillowin zone, except for two symmet- 
rical points where the low energy excitations are gapless 
fermions with a linear dispersion relation^. The low en- 
ergy description of graphene can be done in terms of four- 
component massless Dirac fermions. In this scenario, cu- 
rious behaviors such as quantized conductivitj'^'^i^i^'ii^, 
unconventional integer quantum Hall effect^i^ii^, and the 
Klein paradoxic (unimpeded penetration through high 
and wide potential barriers), can be discussed. 

With regard to conductance, although current exper- 
iments observe a dissipative (Ohmic) regimelS, from a 
theoretical point of view, the study of transport proper- 
ties of charged 2d massless Dirac modes is a fascinating 
problem, far away from its complete understanding^'^. 

In ref. [3, for a thin ultraclean graphene strip, a quan- 
tized conductance related with the zero modes at the 
edge of the strip was found; similarly to what happens in 
nanotubes, where evidence of conductance quantization 
was recently reportedi^. The value of the conductance 
depends on whether the edges are cut with zig zag or 
armchair geometry. In that work, a one particle approx- 
imation was considered, disregarding the effect of any 
kind of impurities or interactions. 

In this work, we claim that this behavior is not only 
correct, but it can be extended to more general situa- 
tions. The conductance quantization is preserved if, in a 
localized region, the system suffers an adiabatic change 
with respect to the situation at the strip, preserving chi- 
ral symmetry. 

We recall that, in planar models involving a four- 
component spinor ip, two independent matrices 7^, 7"* 



that anticommute with the Dirac matrices can be de- 
fined, which play a similar role to the usual 7^, in one 
and three-dimensional Dirac theories. 

At the classical level, a planar system of massless Dirac 
fermions, including (total) charge density and current 
interactions, are symmetric under the continuous chiral 
transformations, generated by 7^, 7*^. Then, as examples 
of adiabatic change, we can consider an adiabatic widen- 
ing of the 2d sample, or even effective current interactions 
coming, for instance, from a smooth deformation of the 
graphene plane^^. 

In the above mentioned conditions, chiral symmetry 
will be sufficient to show that Landauer conductance is 
dominated by the leads, as occurs in Id systems such 
as quantum wires, adiabatically connected to Fermi liq- 
uid reservoirs^^JJ,. Then, following a line of reasoning 
similar to the one presented in refs. ,18. andfiol. a simple 
argument relying on general properties of the system will 
be given, in spite of the fact that the 2d sample contains 
complicated gapless modes, modeled as confined massless 
interacting Dirac fermions. 

From a technical point of view, we will use the func- 
tional bosonization for 2d systems, where the current is 
mapped into a topological current, containing effective 
"electric" and "magnetic" fields for a vector potential 
At, Ax, Ay (the bosonizing fields). In principle, a closed 
form for the dual bosonized action describing 2d Dirac 
fermions is only known in the large mass limit. However, 
on the leads, as the typical scale of one of the dimen- 
sions is small, we argue that results from Id bosonization 
can be applied there. This amounts to approximating, 
on the leads, the 2d fcrmion determinant by considering 
only the contribution of fcrmion zero modes, and taking 
into account the decoupling of the other modes for small 
widths^S. 

Another important ingredient in our derivation will 
be the universal character of the bosonization rule for 
fermion currents in a general 2d system^!. In some sense, 
we will see that bosonization implements the idea of an 
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electron wave guide, enabling the discussion of trans- 
port properties in terms of similar concepts associated 
with "electromagnetic" wave guides, such as geometric 
and material dispersion for the coupled modes. For in- 
stance, the decoupling of the zero modes on the quasi 
one-dimensional leads amounts to the suppression of the 
geometric dispersion, in wave guide language. 

In section i jlll we describe the bosonization technique, 
stressing its physical meaning and showing the similar- 
ities and differences between Id ( Hll Al) and 2d systems 
( mi B[) . In particular, in ijll Bl and mi C[ we explain how 
to construct the bosonized action for the general system 
(2d sample plus quasi one-dimensional leads), showing 
how the bosonization technique implements the idea of 
electron wave guide. In section mill we derive the main 
result of this paper, showing the conductance quantiza- 
tion for a wide class of ultraclean systems with interac- 
tions. Finally, we discuss our results in mVl 

II. CHIRAL TRANSPORT AND 
BOSONIZATION 

A. One dimensional massless fermions 

In this section we give a brief summary of one di- 
mensional bosonization. Although the procedure is very 
well establishedSSi^i^'2^, we would like to emphasize the 
physical concepts involved that will lead us to new re- 
sults, developed in the following sections. 

The description of a conserved charge in a Id system, 
is usually simplified by introducing a bosonic field (p such 
that, 

p = dx(l) , jx = -dt(j), (1) 

which automatically leads to the continuity equation, 

dtp + dxjx = (2) 

for this reason, the 0-field expressions are called a "topo- 
logical current" . 

For instance, we can consider a gapless Id noninter- 
acting fermionic mode ip characterized by a linear dis- 
persion E{p) = vp. In quantum wires, this model can be 
used to represent those modes in the 2c? Fermi liquid in 
the leads that couple to a quasi-one-dimensional quan- 
tum wire (sample), in this case E, and p refer to energy 
and momentum measured with respect to a Fermi surface 
(points), at low excitation energies. In this case the field 
mode can be either right or left moving (ipn — ipn^x — vt) 
OT ipL — ipLix + vt), respectively). Then, the current den- 
sities, 

P = e{PR + Pl) , jx ^ ev{pR- pl), (3) 
where pu = a-ud pL — iplipLi also satisfy, 

v-^ dtjx + dxP = 0, (4) 



in addition to the continuity equation. Then, from eqs. 
(HI) and ([U, the field (p must obey the wave equation, 

{v-'d? ~ dl)cp = 0, (5) 

and we can introduce an action, 

^oM = j dtdx ^ [v-\dt(Pf - [dx^Pf] , (6) 

whose minimization leads to the wave equation ([5]). 

For the associated quantum theories, the quantum 
equivalence between a Id massless Dirac field and a mass- 
less scalar field cj), with the current mapping eq. ([T]) is well 
known and is called bosonization22ii2^'21i2^. For a general 
fermionic theory whose action is Kp[ip] having arbitrary 
short or long ranged interactions /[p, jx] only involving 
densities and currents, we have shown in ref. l2l| the fol- 
lowing equivalence between partition functions, 

Z = j VijjVi) e-y.^iKFW\+iI[p,3x] 

= j'Dcp exp^i^B[0] + il[dx(j). -dtcp], (7) 

where Kb can be computed as a "transverse" functional 
Fourier transformation of the partition function associ- 
ated with Kp (see ref. fiol ). Moreover, the computation of 
Kb for a Id massless Dirac field leads to Xs[(/>] = 5*0 [</>], 
with a = {1/v) In particular, for an external electric 
field, le = J (Px pV, in (/)-language we have a new term 
le = J d'^x{dx(j))V — J d'^x 4>{~dxV), so that the saddle 
point equation reads, 

{v-'d? - dl)cP = aE. (8) 

On the other hand, at the classical level, in addition to 
the usual C/(l) symmetry V'l e'^Vij V'fl: ~^ e'^Vfl lead- 
ing to charge conservation. Id massless Dirac fermions 
with current interactions posses a chiral symmetry, V'l ~* 
e^'^ijjL, 4'R e.~^'^'4'R- Because of the relative minus sign 
in the phases, the right and left mode contributions to 
the chiral density and current are, 

P"^ ^ e{-pB, + Pl) , = -ev{pR + pl), (9) 

containing a change of sign with respect to the U(l) quan- 
tities in eq. that is, 

P^'^-Jx/v , j^^-pv. (10) 

Because of this symmetry, in the classical system we have 
the conservation law, 

dtP^ + d.Jx = 0. (11) 

However, if the action mapping in eq. ([7]) is considered, 
which is needed for the quantum equivalence between 
the fermionic and bosonic partition functions, using the 
(^-language in eq. ([5]), we get, 

dtP^ + d.J^ = viv-^d^ - dl)cp = {va)E. (12) 
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In other words, eq. (jl2p represents the nonconservation 
of the chiral current at the quantum level. A classical 
symmetry that is not realized at the quantum level is 
called an anomaly. From a physical point of view, the 
anomaly represents the creation of particles and holes 
out of the vacuum or Fermi sea, when an external field is 
applied; of course, these concepts only exist in the quan- 
tum world. As it is well known, this effect fixes the value 
of a as follows. Consider an homogeneous electric field 
E and an initial Fermi sea. From the equation of motion 
p = eE all right moving particles will gain a momentum 
eEt. If the system size is L, then the number of right 
moving particles created above the Fermi sea, at pp — 0, 
is the volume occupied in phase space, L{eEt), divided 
by Planck's constant h, so that the density of right mov- 
ing particles is eEt/h. This equals the density of left 
moving holes (or antiparticles) created, to conserve the 
system's charge equal to zero (with respect to the charge 
of the Fermi sea). This is measured by the quantum ver- 
sion of p = p-{pb. + Pl)- Now, because of the relative 
minus sign in pA = e{—pj^ + p^), at the quantum level, 
it will measure the total density of particles (times the 
electric charge e): pA = {^e^Et/h). Since the system is 
homogeneous, eq. (2) then gives, 

dtPA + d^j^ - 2e'ElK (13) 

and the value, 

a = {l/v) ^ (14) 

is obtained. Note that if N dispersionless one- 
dimensional channels ipi, i = 1,..,N were considered, 
each one represented by a field 4'i, the anomaly for the 
total chiral currents in eq. would be {vNa)E, a = 
{l/v) Equivalently, if a single field were used to de- 
scribe the total currents, then the value a = (l/v) 
should be used for the associated field theory. 

B. Two dimensional four-component massless 
fermions confined to a strip 

Following the same reasoning of the previous section, 
let as consider a 2d system with a conserved charge, and 
the associated continuity equation, 

dtp + dj;jx + dyjy = 0. (l5) 

Formally, this equation is a "three-divergence" equal to 
zero, so that we can represent p, jx and jy as a "curl" or 
topological current, 

p — dxAy — dyAx 
jx = dyAt- dtAy 
jy = dfAx - dxAt, 

(16) 



which is identically conserved. 

Considering that At, Ax, Ay and At + dtX, Ax + dxX, 
Ay+dyX represent the same charge density and current 
distribution, the A-fields must be gauge fields, physically 
equivalent when the above mentioned transformation is 
performed. It is useful to think about the charge density 
and currents as effective "magnetic" and "electric fields" . 
Making the identifications, we have, 

p ^ B 

jx = +Ey 

jy = ^Ex- 

(17) 

Of course, E and B are not real electromagnetic fields, 
they are simply useful auxiliary fields to represent charge 
density and currents. In particular, charge conservation 
now looks like a 2d Faraday-Lenz law, 

dxEy - dyEx - -dtB. (18) 

Thus, charge variation in a given region is associated with 
nonzero current flux through the boundary, as changes in 
the magnetic flux piercing a surface are associated with 
an induced electric field, in bosonized language. 

If gapless parity preserving fermions 'ijj with a linear 
dispersion E{p) = v\p\ (no "material" dispersion) are 
considered, we can guess that in 2c? we will have to face 
a difficult problem because of "geometrical" dispersion. 
The fermion modes are general combinations of waves 
eicpi{k.x ± Lot) propagating with speed v — uj/\k\, satis- 
fying the wave equation derived from Dirac's equation. 
But the charge density and currents are i/j-field bilinears, 
not satisfying any simple equation, because of the con- 
tinuum of possible directions given by k ("geometrical" 
dispersion). That is, if we try to write other equations 
defining p, and j, or their effective electric and magnetic 
counterparts, they will certainly be highly nontrivial. 

On the other hand, if gapless excitations were confined 
to a 2c? strip, we expect that the effects of geometrical dis- 
persion will decrease, as the strip width is reduced. Then, 
we could use the (/)-language representation of miA\ for 
these modes. However, we would like to consider a gen- 
eral situation where the strip could be coupled adiabati- 
cally with an extended 2c? region, where a continuum of 
gapless modes exist, with no definite direction of propa- 
gation, described by the gauge field A. For this reason, 
we will translate the dispersionless modes in the strip 
from (j) to ^-language. 

If the strip is defined by a region of length L along 
the X-axis, and width W, that is y e [~W/2, +W/2], the 
confinement condition of having no charge flux across the 
limits, that is jy = at y = ±W/2, amounts to Ex = 
at y — ±W/2. In terms of the effective model this means 
that the the border acts as a perfect conductor or wave 
guide (see ref. [26h . 

Assuming that the current distribution on the strip is 
all along the x-axis, that is jy = 0, while p and jx are 
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y-independent, we have, 

dxjy - dyjx = , w"^ dtjy + dyp = 0, (19) 

or in "effective" ianguage, 

dxEx+dyEy = , dyB~v-^dtEx^O. (20) 

Now, using tlie condition ([4]), for dispersionless Id modes, 
we also liave, 

dxB + v-^ dtEy = 0. (21) 

Eqs. mi), dlOl) and are Maxwell's equations in 2d, 
while the current distribution considered in the strip cor- 
responds to an effective dispersionless TEq mode coupled 
in a wave guide. These equations can be derived from 2d 
Maxwell's action, 

So[A] - J dtd^x^[v~'{El + El) - B% (22) 

with Ex, Ey and B defined through eqs. (|16p and P7)l . 

We also note that in the radiation gauge dxAj;+dyAy = 
0, At = 0, and for the TEq mode we can use Ax = 0, so 
that eq. p6|) reads, 

p = +dxAy 
jx = -dtAy. 
Jy = 0. 

(23) 

Then, considering that p in Id and 2d correspond to 
charge by unit length and area, respectively, while jx 
has units of current and current by transverse length, 
respectively, we can identify Wp and Wjx in 2d with 
one dimensional quantities, that is, we can identify (j) = 
WAy. We would like to underline that the dispersion- 
less modes have a formal "Lorentz" invariance associated 
with boosts with speed parameter v. Under this boost Ay 
is invariant, so that, in the quasi one-dimensional lead, 
its identification with the scalar </> is a natural one. 

Summarizing, if we evaluate Maxwell's action on the 
TEq mode, we obtain the identification /3 — a/W, 
a — (l/v) for the description of the total currents 

associated with N dispersionless channels on a 2d strip 
be equivalent to N Id dispersionless quantum modes. 

C. General chiral symmetric fermions confined to a 
2d sheet 

For a general fermionic 2d model with interactions 
I[p, jx , jy] , only involving densities and currents, the gen- 
eral quantum equivalence with the ^-language descrip- 
tion is given by, 

Z = y VipV-tp exp iKp [iIj] + il[p, jx , jy] 

= J VAexpiKB[A]+iI[B,Ey,-Ex], (24) 



where Kpltp] is the free fermionic action. The bosonizing 
action is gauge invariant, and again is given by a 

"transverse" Fourier functional transform of the partition 
function associated with Kp^^. 

The low energy description of graphene can be written 
in terms of a four component spinor ip, which can be or- 
ganized as a pair ipx, V'-ifi where the label is associated 
with excitations in the two opposite valleys centered at 
the corners of the Brillouin zone, with wave vector ±K 
(Dirac points). The fields 'iI'±k are two-component Dirac 
spinors, each component is associated with the amplitude 
of the wave function on the A and B sublattices of the 
honeycomb lattice. 

We underline that while the effective action for two- 
component massless fermions contains a Chern-Simons 
term^^, in the case of four component fermions, each ele- 
ment of the pair must contribute with an opposite sign, in 
order to conserve parity. Then, as the model for graphene 
is based on a four-component ijj, Kb contains no Chern- 
Simons term, and being gauge invariant it must be nec- 
essarily of the form, Kb = Kb[B, Ex, Ey]. 

In general, if Kp were a gapped (massive) theory, it 
would be possible to obtain a low energy expansion for 
Kb- However, in graphene the excitations are associ- 
ated with massless Dirac fermions, and in an extended 
2d region is an unknown complicated functional, 

as there is no parameter to organize a perturbative ex- 
pansion. 

On the other hand, from the previous discussion, on 
a thin 2d strip, iirB[v4] is expected to have the simple 
Maxwell form of eq. (|^^ . as the problem can be con- 
sidered as quasi one-dimensional. To be more precise, a 
massless theory confined to a thin strip will have zero 
modes plus other modes that decouple in the small W 
limit (see ref. |20|). Thus, the partition function for Kp 
on the strip is expected to be dominated by these zero 
modes; these are precisely the relevant modes to discuss 
transport in ultraclean graphene strips (see ref. 13). 

Then, in the more general case, where a two di- 
mensional system is adiabatically coupled to quasi one- 
dimensional leads, the mapping (|24p can be used, if cor- 
rections to the simple Maxwell form in the leads are in- 
corporated in the bulk expression for Kb- 

III. CONDUCTANCE OF A GENERAL 
GRAPHENE STRIP/SAMPLE/STRIP SYSTEM 

From the previous discussion, if an extended 2d region 
with ultraclean graphene (gapless Dirac fermions) is adi- 
abatically coupled to two quasi one-dimensional strips, 
formed with the same material, the system's bosonized 
action is expected to have the form, 

Sb^Kb + I , Kb = Sq + R, (25) 

where 5*0 is given by eq. ()22p . R represents deviations 
from the simple Maxwell form valid on the leads, local- 
ized in the sample's bulk. The /-term represents charge 
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density and current interactions also localized at the bulk 
In the present bosonized language, these interactions and 
the deviations R play a similar role. In addition, the 
"electric" field E^, Ey must be orthogonal to the sam- 
ple's boundary^^ (no charge flux across the boundaries). 

The important point is that because of the adia- 
batic condition, Kb is defined by a single functional 
of B, Ex and Ey, whose form depends on the system's 
(lead/sample/lead) point under consideration. Although 
the corrections in the bulk are very difhcult to compute, 
we will see that the above mentioned general structure 
for Kb is all we need to derive the conductance quanti- 
zation. 

To discuss transport, we include in eq. (j24p a sim- 
ple coupling with external electric and magnetic fields 
{Ao,Ax,Ay} to probe the system, that is, we consider 
the replacement, 

I ^I + Ie , Ie= J d^X {~pAo + jxAx + jyAy), 

(26) 



or in A-language, 



1^1 + j dPx {EyAx ~ ExAy - BAo 



(27) 



The associated saddle point equations are obtained by 
taking functional derivatives with respect to At, Ax and 
Ay, respectively, 

8Sb , „ 5Sb „ 



SEx 



5E, 



n 6Sb , ^ 5Sb c- 



SB 



6Ex 



5Sb SSb r. 



SB 



6E„ 



(28) 



where we have defined the system's bosonized action, 



Sb=Kb + I. 



(29) 



Now, we can follow a reasoning very similar to the 
one of Maslov and Stoneil used to derive the universal 
behavior of conductance in quantum wires (see also ref. 

Let us consider an electric field that is switched on 
to attain a stationary value £x{x,y), Ey{x,y). At late 
times a uniform current / is expected to be settled in the 
system, so that on the strips we can consider an ansatz 



0, Ax 



and Ay 



-dtAy — k, so that / = kW. 



fix) - kt, I/W = jx 



Because of causality. Ay on the left (right) lead must 
be given by a perturbation propagating to the left (right), 
with speed v. Then, considering symmetry under x- 
reflection (to simplify the argument), /(x) must have the 
form —{k/v)x on the left lead and +(k/v)x on the right 
lead. This ansatz can be extended to a solution on the 
whole system, and considering the last two equations in 



28]), we have. 



dx 



-(1//3)S + 



5{R + I) 

SB 
S{R + 1) 



6B 



= -£y 



(30) 



Integrating on a curve C going from the left to the right 
lead, 

- {l/(3)[B\„ght - B\ieft] = ~ j [dxSx + dySy) = A^, 

(31) 

and using eqs. and (fT7|) . and the ansatz above. 



B\right — {,^xAy dy Ax) right — 



B\left — {dxAy — dyAx)left 

Therefore, we obtain, 



(32) 



/ = ~(Wv(3/2)AV = -N—AV, (33) 

where N is the number of zero fermion modes on the 
strips. 



IV. SUMMARY AND DISCUSSION 

We have analyzed the problem of linear transport in 
an arbitrarily shaped ultraclean graphene sheet, adiabat- 
ically connected to quasi one-dimensional leads, formed 
with the same material. This means that a smooth 
widening of the strips is permitted at the sample's bulk, 
where effective interactions are localized, which depend 
on the total charge density and currents associated with 
both Dirac points. Under this circumstances, we have 
shown that the conductance is quantized, as it is domi- 
nated by the ultraclean quasi one-dimensional leads. 

The situation is similar to what happens in quantum 
wires, where conductance is quantized because of the 
dominance of the Fermi liquid reservoirs. 

There is an interesting realization of the scenario dis- 
cussed here, recently studied in the literature, quite rele- 
vant when one tries to compare models with actual exper- 
iments. In ref. [l^ Castro Neto and Kim analyzed the ef- 
fect of the curvature of the (otherwise plane) carbon sheet 
on the electronic structure of the sample. They showed 
that the net effect of smooth curvature fluctuations is 
that fermions become minimally coupled with an effec- 
tive "elastic" gauge field. Permanent deformations could 
eventually break time reversal symmetry at one point. 
However the effect on the other Dirac point —K restores 
this symmetry. In our formalism, that means that the 
bosonized action Kb has no induced Chern-Simons term, 
due to the cancellation between both species of fermions. 
On the other hand, the elastic field couples symmetrically 
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the two Dirac points, meaning that smooth fluctuations 
of this gauge field will induce total charge and current 
interactions of the class considered in this paper, so that 
our quantization argument can be applied in this case. 

In general, our argument works whenever the chi- 
ral symmetry displayed by 2d four-component massless 
fermions is preserved, so that the obtained perfect con- 
ductance can be associated with effective gapless exci- 
tations defined along the whole system. Bosonization 
gives a nice interpretation for this phenomenon, looking 
to the sample as an "electromagnetic" wave guide, with- 



out backscattering, guiding the TEq modes at the leads. 
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